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We study a composite quantum quench of the energy gap and the interactions in the interacting 
0^ model using a self-consistent approximation. Firstly we review results for free theories where a 
quantum quench of the energy gap or mass leads for long times to stationary behaviour with thermal 
characteristics. An exception to this rule is the 2d case with zero mass after the quench. In the 
composite quench however we find that the effect of the interactions in our approximation is simply 
to effectively change the value of the mass. This means on the one hand that the interacting model 
also exhibits the same stationary behaviour and on the other hand that this is now true even for 
the massless 2d case. 
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I. INTRODUCTION 

An area that has been gaining increasing interest over 
the last years is that of out-of-equilibrium quantum 
physics. An example of particular simplicity is that of 
quantum quenches in which some of the parameters of the 
hamiltonian of an isolated quantum system are changed 
instantaneously. Then one practically has to study the 
time evolution of a trial wavcfunction, which is typically 
the ground state of the hamiltonian before the quench, 
under the influence of the hamiltonian after the quench. 
Although one expects a periodic collapse and revival of 
the initial state, in practice this period diverges rapidly 
with the system size and for large systems local observ- 
ables may exhibit stationary behaviour at long times, 
eventhough the global wavcfunction itself may never be- 
come such. This has been shown to be the case in many 
different settings^"— . 

An obvious interesting question is whether this station- 
ary behaviour is thermal as one may reasonably expect. 
It turns out that in many intcgrable systems the sta- 
tionary behaviour is described by a statistical distribu- 
tion which is similar to but not exactly thermal^"— i^ d^i^^ . 
More precisely it is a generalized Gibbs ensemble, subject 
to the constraints imposed by the integrals of motion. It 
was then conjectured that non-integrability is responsi- 
ble for the exact thermalization of a systeroii. To what 
extend this is true is however still under investigation, 
since theoretical arguments that support this conjecture 
are based on semiclassical conjectures^^, while numeri- 
cal studie a^^d'''^^' '— lead to rather controversial results: 
some of theroi^ii^i2£ reveal non-thermal behaviour even 
for non-integrable systems, while othersiii^ are in good 
agreement with the thermal predictions and attribute 
the previous disagreement to finite size effects. There 
are some analytical studies in lattice models too^^: the 
firsb^ refers to the Bose-Hubbard model but after the 
quench the system evolves under the free hamiltonian 



of the superfluid regime. In the secondS^ dynamical 
mean-field theory (DMFT) applied to the non-integrable 
Falicov-Kimball model shows non-thermal features. On 
the other hand, an interaction quench in the Fermi- 
Hubbard model is possible to lead to thermalization, 
as shown using two different analytical^! and numerical 
DMFTi2^ approximations. 

In the present work we will study quantum quenches 
employing a field theoretic approach, which is supposed 
to capture their essential general characteristics. We con- 
sider systems described by a relativistic dispersion rela- 
tion with some energy gap (or mass) and a maximum 
group velocity of excitations. Then for free systems, a 
quantum quench of the energy gap leads to stationary 
behaviour and a momentum dependent effective temper- 
ature can be defined^i^. This is true for quite general con- 
ditions: the energy gap after the quench must be nonzero 
in Id and 2d while in 3c? the result holds even if it is zero. 
Furthermore a, Id gapless system can only be interacting 
and it turns out that it exhibits similar behaviour too^i^. 
It should be emphasized that the notion of effective ther- 
malization used throughout the present and some related 
earlier work^ refers to the thermal-like stationary be- 
haviour that fits to the generalized Gibbs ensemble de- 
scription rather than the standard thermal theory. This 
is manifest in the fact that each momentum mode corre- 
sponds to a different effective temperature, since in the 
absence of interactions each mode evolves independently 
from the others. This suggests that in interacting systems 
the energy exchange due to collisions between different 
momentum modes would result in a mixing of their effec- 
tive temperatures. However if the interaction is such that 
the system is still integrable then there will be some other 
decomposition into independent modes (quasiparticles) 
and we expect that the system still exhibits stationary 
behaviour with a different effective temperature for each 
of these modes. Therefore it is only when the interaction 
makes the system non-integrable that thermalization to 
a unique common temperature is still a possibility. 
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The simplest interacting field theory is the (j) model. 
We consider a simultaneous quench of the mass from tuq 
to m and of the coupling constant from Aq to A. In 
order to study the evolution of the system we need to 
use an approximation scheme and the simplest one is 
the Hartree-Fock or self-consistent approximation. This 
can be applied in a number of different but equivalent 
ways. In perturbation theory it consists in ignoring all 
skeleton diagrams from which the diagrammatic expan- 
sions of correlation functions are constructed, except for 
the simplest one, i.e. the loop diagram. This turns out 
to be the same as approximating the system's state by 
gaussian wavefunctions or substituting the quartic in- 
teraction term in the hamiltonian by a quadratic one 
with a self-consistent coefficient. Notice however that in 
this simple approximation, collisions between particles 
of different momenta are neglected and this makes our 
approach incapable of answering the previous question 
about the relation between non-integrability and exact 
thermalization. Indeed, although the t/)^ model is non- 
integrable, the Hartree-Fock approximation becomes ex- 
act only in the large- A'^ limit of the linear cr-model, i.e. 
the generalization of the 0^ model to an A'^-componcnt 
field, which becomes integrable in this limit. Thus our 
approach provides the integrable counterpart of a non- 
integrable model that best approximates it. It is however 
the necessary first step towards understanding the effect 
of quantum quenches in interacting systems and should 
be expected to reveal some of their general features. 

There is a significant number of publications that use 
the same method to study other closely related out-of- 
equilibrium problems, partially due to applications to 
cosmology. Cooper and Mottola3^ make a detailed pre- 
sentation of the method for the evolution of a general trial 
wavefunction and Boyanovsky ct ali^"— study the spe- 
cial case of a quench from the disordered to the ordered 
phase at large temperature and in 3d. Also Wetterich et 
al. have studied the time evolution of out-of-equilibrium 
initial ensembles using a different method based on the 
numerical computation of the time-dependent effective 
action^Si^. Recently a remarkable numerical study based 
on the same method and including next-to-leading order 
effects in the large- iV expansion, has shown that an ini- 
tial pure state evolves so that the reduced density matrix 
indeed thermalizes at large times^. 

Using our approximation we find that the two point 
correlation function long after the quench is of the same 
form as the free correlation function but with a different 
mass that has to be determined self-consistently. This 
means that nothing really changes in terms of the relax- 
ation of the system: once again it becomes stationary and 
a momentum dependent effective temperature can be de- 
fined, the only difference being that m will be replaced by 
an effective mass m* which depends also on the coupling 
constant A. The self-consistency equation for m* has al- 
ways a real solution larger or equal to m. In the critical 
case m = 0, we find that in Id m* is also zero, but in 2c? 
it becomes finite. This leads to the important conclusion 



that in 2d, after a quench to zero mass which according to 
the above discussion would not lead to relaxation if the 
system were free, now due to the presence of the inter- 
action, it acquires a non-zero effective mass which allows 
it to relax. Furthermore, by studying the time evolution 
of the effective mass wc find that if mg > m and A is 
sufficiently large then right after the quench the system 
is effectively set into an unstable state although it soon 
recovers its stability. 

In the first part of this paper we focus on free sys- 
tems which have been partially discussed earlier— i^. Here 
we present an elegant simplified derivation of the quench 
propagator, develop an exact imaginary time formulation 
based on an earlier invented mapping to a slab geometry 
and define an average measure of the effective tempera- 
ture first introduced in recent work2^. These constitute a 
useful toolkit for many applications and extensions. For 
completeness we briefly report earlier results regarding Id 
integrable systems with critical evolution. In the second 
part we study the composite quench in the model in 
the self-consistent approximation. This part is split into 
two sections: in the first we follow a heuristic approach 
based on perturbation theory and find an ansatz for the 
correlation function and in the second we start with the 
equations of motion and investigate the time evolution 
to verify the results obtained from our ansatz. 

II. SIMPLE HARMONIC OSCILLATOR 

The simplest problem of a quantum quench one can 
start with is that of a simple harmonic oscillator whose 
frequency is quenched from wq to w. The hamiltonian 
before the quench is 

Ho = \-K^ + (1) 

while after the quench it is 

H=\.' + \y<i>' (2) 

The initial state is the ground state |^'o) of Hq. 

From a physical point of view, what happens is that 
I'I'o), as a trial state different from the ground state 
|0) of H , contains, compared to that, an energy excess 
which is distributed to the excitation levels of H . Af- 
ter the quench the evolution of the wavefunction in the 
Schrodinger picture is given by 

|*(t)) = e-*^*|*o) = ^e-'("+i/2)"*|n)(n|*o> (3) 

where \n) is an arbitrary eigenstate of H . 

It is trivial to observe that the evolution is periodic 
since after a period T = 2tt/uj the system returns back 
to the initial state, up to an irrelevant minus sign. This is 
a special case of quantum recurrence^. In fact the wave- 
function will exhibit periodicity or quasi-periodicity (i.e. 
it will return arbitrarily close to the initial state after 
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sufficiently large time) in any system with discrete en- 
ergy eigenvalues. Systems with finite degrees of freedom 
always have such discrete spectra, while in the thermody- 
namic limit the spectrum becomes in general continuous 
and quantum recurrence may be lost. In practice even 
for finite but large systems, the corresponding period is 
usually so large that this periodicity is irrelevant. 

A. Propagator 

We are also interested in the correlation function of 
the field operator at different times, i.e. the propaga- 
tor (*o|7'{^(ti)<?!)(i2)}|*o) = Cq{ti,t2) where T denotes 
time ordering. The time evolution of x in the Heisenberg 
picture is given by the equations of motion 

+ = (4) 

which can be solved easily 

, , , , , sin Lot , , 

= 0(O)cosa;t-t-7r(O) (5) 

We therefore have 



Notice that we have separated the Feynman propagator 
^-iui\ti-t2\ which, as expected, is the only term that 
survives \i uj ~ bjQ i.e. if there is no quench at all. Also 
notice that the only term that breaks time invariance is 
the second one. 



III. LINEARLY COUPLED OSCILLATORS 
(FREE FIELDS) 

Let us now move on to study a system of linearly cou- 
pled harmonic oscillators or equivalently a free field the- 
ory. In general such a system is described by a quadratic 
hamiltonian of the form 

ii-\Y.^\r) + \Y^K(r- r'^r) ^(r')f (9) 

which can be easily diagonalised in momentum space 
where it takes the form 

^ = Y^ ^^fe^-fc + \^l'i>k^-k (10) 



(Vl/o|0(ti)0(t2)|*o) 

+ (*ok'(0)|*o) 



*ol0^(O)|5'o) COSwti COSLOt-i- 

sinwti sinwt2 



(*o|'/'(0)40)+ 71(0)0(0)1*0 



sina;(ti + ^2) 

. sina;(ti - t2) 
* 2aJ 



(6) 



where the canonical commutation relation [0(0), 7r(0)] = 
i has been taken into account in order to simplify the 
last term. All terms are symmetric under the interchange 
t\ O ti apart from the last one which is antisymmetric. 
Thus time ordering amounts to substituting {t\ — ti) in 
the last term by its absolute value. 

It is now clear that the problem reduces to the calcula- 
tion of the initial expectation values of 0^(0), 7r^(0) and 
(/)(0)7r(0) -I- 7r(0)(/>(0). From the initial condition that the 
system lies in the ground state of iJo we easily find 
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(*o|0(O)7r(O)+^(O)0(O)|«'o) 



(vI/o|02(O)|vI'o) = ^ 
(*ok'(0)|*o) = 



= 



(7a) 

(7b) 
(7c) 



and substituting into © we obtain 
Cq{tiM) = ^^,^2^!^°^ cosu;(ti-^2)+ 



^cosw(ti +i2) + 



_}_^-iui\ti-t2\ 
2UJ 



(8) 



We will assume a relativistic dispersion relation 
2 2 7 2 I 2 4 



(11) 



with energy gap (or mass, in the language of quantum 
field theory) m and speed of sound c. This can also de- 
scribe successfully non-relativistic systems with the same 
energy gap m and maximum velocity of excitations c. 

The quantum quench that we will consider consists in 
an instantaneous change of the mass from mo to m. For 
brevity we can set c = 1. An investigation of a quench 
of the speed of sound c is done elsewhere^. As earlier, 
we assume that before the quench at t = the system 
lies in the ground state of the initial hamiltonian |^'o). In 
addition the system is kept isolated from the environment 
before and after the quench. 

In order to study the time evolution, it is sufficient to 
find the two-point correlation function i.e. the propaga- 
tor 

(*o|r{0(ri,ti)0(r2,i2)}|*o) =C,(ii,t2,ri-r2) (12) 

since in a free theory all physical obscrvables can be ob- 
tained from this. From (jlOp we see that the system is 
decomposed into a set of independent momentum modes 
each of which evolves as a simple harmonic oscillator. 
Thus the propagator is simply the Fourier transform with 



respect to k of the expression 
and (jJk = \/k^ + -m? 



Cq{tl,t2,r) 



cfik 
{2tj:Y 



with cjofe = -^/fc^ + ; 



e*''C,(ti,t2;fc) (13) 
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A. Properties of the propagator 

Let us now study the physical properties of the equal 
time propagator in real space. For simplicity we will 
mainly use its asymptotic form for mg 3> m and t, ?■ 3> 
tUq^- This will be called the deep quench limit and should 
obviously exhibit all characteristic features of a quantum 
quench since it is one of the two most extreme possibili- 
ties for the relation between the two masses. In this limit 
the propagator simplifies to 



d'^k Too 



e""' 7-3(1- cos 2L^fei) (14) 



The masslcss (to = 0) and massive (to 7^ 0) cases are 
different and should be investigated separately. 



of the integrand we close the integration contour in the 
upper half plane but for the time dependent part we have 
to rotate it by 90° instead. Each part has single poles 
at fc = ±im and the outcome is nontrivial. Exact re- 
sults cannot be found and we have to employ asymptotic 
methods for large r and t. In particular using the sta- 
tionary phase method we find that for fixed r and large 
t the time dependent part of the integral tends to zero 
like t"''/^ cos 2mt. Also the rest decreases for large r like 

g-mr^^(d-l)/2^ 

We thus conclude that the propagator changes sharply 
as we cross the lines t = r/2. Before this time there are no 
correlations between two distant points, while afterwards 
the two points become correlated. This feature, which is 
a direct consequence of the causality principle, is called 
the horizon effect. Fig. [1] illustrates the main features of 
the massive propagator in Id. 



1. Massless case 

In this case = and after some algebra using 
Fourier transforms of common functions, we obtain the 
following exact results: 



d=l 




r)/8 



if r > 2t, 
if r < 2t. 



(15) 



• d = 2 

C!tf\r,t) 



!^\og[{2t + VW^)/r] iir<2t 



if r > 2t, 
It. 
(16) 



d = 3 







if r > 2t, 



^S'^(^'^)-C^/16^,. |fri2Z ^^^^ 



In all dimensions we distinguish between two spacetime 
regions in which the behaviour of the propagator is qual- 
itatively different: for r > 2t it is always zero, unlike 
for r < 2t. This means that the correlations between 
two points at distance r remain unchanged until t = r/2. 
Also notice that in id the propagator is time independent 
for r < 2t. 



2. Massive case 




C{ro,t)/mo 




FIG. 1: Top: Spacetime plot of the deep quench propaga- 
tor Cdq{r,t) in Id and for m = 1, as obtained by numeri- 
cal integration of (|14|l . The horizon effect is clearly demon- 
strated. Outside the horizon the value is exactly zero. Bot- 
tom: Time dependence of Cdq{r,t) (blue hue) at fixed dis- 
tance r = ro = 2, denoted by the vertical red line in the 
above figure. The dashed lines give the large time asymptotic 
expressions. Notice the decaying oscillations ~ t~^^'^ cos2mt 
(purple line) around the stationary value ~ e"™*^ (red line). 



By evaluating the integral ([H)) we notice that, as be- 
fore, we have to distinguish between two spacetime re- 
gions in which the behaviour of the propagator is qual- 
itatively different. If r > 2t then we can close the inte- 
gration contour in the upper half of the complex fc-plane 
and since there is no pole the integral is zero. If on the 
other hand r < 2t then for the time independent part 



Another particularly important conclusion is that if 
m then for fixed distance the propagator becomes 
stationary for large times. The same is true for m = in 
3d, but not in Id or 2c?. In addition this result is robust 
and does not rely on the deep quench approximation. 
Indeed if we use the full expression of Cq{t;k) ^ for 
m — and 3d we find that the time dependence decays 
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exponentially. The Id case is more complex and requires 
special treatment. We will talk about this in section HVl 



B. Comparison with the slab propagator 

We will now study a completely different problem 
which however turns out to be an imaginary time for- 
mulation of a quantum quench. We consider a euclidean 
free field theory defined on a (c? + l)-dimensional slab 
of thickness L with Dirichlct boundary conditions, that 
is the two-point correlation function or Green's function 
vanishes when one of the points are on the boundaries of 
the slab r = — i/2 and t = +L/2, where r is the trans- 
verse coordinate. The Green's function G'si(ri, ri, r2, T2) 
for this problem can be found using the method of im- 
ages as follows: to reproduce the boundary conditions we 
put an infinite set of alternating positive and negative 
'charges' at the reflections of the 'source' on the bound- 
aries (Fig. Oa). Then Gs;(ri, ti, r2, T2) is the superposi- 



(a) 



(b) 



By analytically continuing to real times r — ^ it we find 



_ cosiUfc(ti - ^2) e"''^ cosa;fc(ti + ^2) 1_ 

" Wfe(e2-fci - 1) " Wfc(e2"'^i - 1) 2wfc ' 



-iujk\ti—ti\ 



(20) 



If we now compare the slab propagator in real time ([20)) 
with the quench propagator (|8]) we notice that these are 
exactly equal if and only if 



1 



(21a) 
(21b) 



Remarkably, the above two conditions are consistent and 
the solution is 



tanh (ajfei/2) 



^ak/^k 



if Wfe > cjofc- 



(22) 



FIG. 2: Images required for the slab with Diriclilet (a) or 
periodic (b) boundary conditions. 

tion of (euclidean) Feynman propagators between (r2, T2) 
and each of the images of the source (ri,Ti). Since the 
problem is translationally invariant in the d longitudinal 
directions, in the mixed (A;, r) representation we find that 
Gsi{Ti,T2;k) is 



\n— 71—1 



—uJk{~\Ti—T2\+2nL) 



T1=0 



-Wfc(Ti+r2 + (2n+l)L) 



^ g-a;fc(-ri-T2 + (2n-l)L) j 
ri=l / 

(18) 



This is a geometric series and the result is 

g-<^fc|ri-r2| ^ g+^fc(|ri-r2|-2L) _ 26^'^'''^ COshLJfc(Ti + T2) 



2cjfe(l - e-2"loi) 



coshwfc(Ti — T2) e'^'' coshci;fe(Ti + T2 



(19) 



Notice that if we solve with respect to L, the answer is a 
function of k. 

Thus the problem of a quantum quench can be equiv- 
alently formulated as a euclidean theory on a slab with 
momentum dependent thickness. The initial conditions 
in real time are translated into boundary conditions on 
the slab. In the deep quench limit mo — >■ 00 the con- 
dition becomes L 2/mo, independent of k and there- 
fore the analogy between the quantum quench and the 
slab is asymptotically exact. The reason is that Dirich- 
let boundary conditions correspond to vanishing initial 
value of the quench propagator, which is indeed the case 
for TTiQ 00, since Cg(0,0; k) ~ l/2ajoA: 0. 

It should be mentioned that our choice of Dirichlet 
boundary conditions has nothing special: in fact it is 
only important in the deep quench limit. One can ver- 
ify that the quench propagator can be similarly identi- 
fied with the Green's function corresponding to the fol- 
lowing general boundary conditions (known as Robin or 
'impedance' boundary conditions due to their application 
to electromagnetics) 

n ( i.\ , ,dGsi{Ti,T2]k) 

aGsi[Ti,T2;k) + b = (23) 

on 



where 

a^ujh smh{uJk/ojQk) - ^ofe cosh(wfe/a;ofe) 

b = cosh(wfe/a;ofe) - i^ok/^^k sml-i{uJk/i^Qk) (24) 

d/dn denotes the normal derivative at the boundary r = 
±L/2 and in this case L is chosen to be L = 2/a;ofe- 
Note that for wofe » uJk the latter condition reduces to 
Dirichlet type. To intuitively understand the meaning of 
these boundary conditions, we can use an analogy from 
electromagnetics. There, Dirichlet boundary conditions 
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correspond to complete reflection by a perfect conductor, 
while Robin boundary conditions correspond to partial 
reflection and refraction by an imperfect conductor with 
a large refractive index. 

The correspondence between a quantum quench and 
the slab construction turns out to be valid, at least in 
the deep quench limit, even in interacting theories where 
an exact solution may not be possible^. 



C. Comparison with the thermal propagator 

Let us now compare the above two propagators with 
the thermal or Matsubara propagator, which describes 
a system at thermal equilibrium at finite (inverse) tem- 
perature /3. As is well-known, in imaginary time this 
corresponds to the Green's function in the geometry of 
a {d + l)-dimensional cylinder of circumference /3, i.e. a 
slab of equal thickness with periodic instead of Dirichlet 
boundary conditions. Among other ways, this can also be 
derived using the method of images. To reproduce the 
periodic boundary conditions we now need to put only 
the positive images (Fig. [21b) and the result is 



(25) 

or in real time, after the analytical continuation t — > it 

(26) 

We observe that if we could ignore the (ti + t2)- 
dependent part then the slab propagator Gsi{ti,t2; k) 
and the quench propagator Cq{ti,t2]k) would be the 
same as the thermal propagator Gth{ii,i2]k) with L ~ 
13/2. This can actually be correct for the real space form 
of the quench propagator at large times, as we have al- 
ready seen in section IIII Al Indeed this is the case if 
m 7^ or m = and d = 3. 

As a conclusion, at large times the system tends to 
a state with thermal-like correlation functions, which is 
what we named effective thermalization. The effective 
temperature is given, according to all the above, by the 
condition 



tanh(/3cffa;fc/4) = 



iOk/ujok if Wfe < wofc, 
ujok/^k if > LOok- 



(27) 



Notice that the effective temperature is momentum de- 
pendent, which could be expected since, as we already 
mentioned, in a free system each momentum mode 
evolves independently from the others and there is no 
reason why they should all thermalize to the same tem- 
perature. Yet in the deep quench limit the effective tem- 
perature becomes momentum independent /3off ~ 4/mo. 

It should be emphasized that the state itself is neither 
thermal nor stationary: the density operator still exhibits 



oscillating behaviour for example. However, since in a 
free system all local observables can be derived from the 
two-point correlation function which does become sta- 
tionary, the same happens to all such observables as well. 
It is crucial that the system is in the thermodynamic limit 
and the observables under consideration are local since 
then an integration over an infinite set of momenta is 
required and it is exactly this interference of all inde- 
pendent momentum modes that leads to thermalization. 
Such observables include those defined on any finite sub- 
system A of the whole system, like the reduced density 
operator of A^. In this sense the complement of A acts 
as a thermal bath with which A comes into thermal equi- 
librium. This explains why the effective thermalization 
that we consider docs not contradict with the fact that 
in a free or more generally integrable system, there is an 
infinite set of conserved quantities that prevent the sys- 
tem from thermalizing as a whole. The subsystem A is 
not closed and there are no such restrictions to prevent 
its thermalization. 



D. Estimation of the effective temperature from 
the field fluctuations 



As we saw, the effective temperature in our free model 
is different for each momentum mode. Since the low mo- 
mentum modes are those that determine the large dis- 
tance behaviour, for most purposes Pcttik = 0) is suffi- 
cient in order to macroscopically describe the system. We 
can define^ however an estimate of the effective temper- 
ature that averages over all momentum modes in a nat- 
ural way, by comparing the field fluctuations long after 
the quench {(l>^{x = 0, t -> oo)) with those of a system at 
thermal equilibrium. We can call this average effective 
temperature and denote it as /3. Then /3 must satisfy 



d'^kC:{k;m,mo) ^ d'^k Gth{k;mJ) (28) 



where C* stands for the stationary part of the quench 
propagator. More explicitly 

^k'^-'dk ^'^^"-'"!:^' = Hk^^-'dk ^ 



(29) 

from which we can find /3 as a function of m and mo. 
The latter can be written in dimensionless form as 



fd{m/mo) = m gdWm) 



(30) 



where 



and 



9d{s) = I k^-^dk 



yp-7l(e«Vfc^ - 1) 



(32) 
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In units of uiq, setting x = m/mQ and y ~ /3mo wc have 



fdjx) 
9d{xy) 



(33) 



Tables U and HIl show the asymptotic behaviour of the 
integrals in several limits for the relation between the 
parameters and when possible their exact form. 

Figure [3] shows a plot of /3 as a function of m in units 
of Too as obtained numerically from the above equation. 
Note that for m — mo i.e. no quench at all, the effective 
temperature 1//3 is zero as it should be. Apparently in 
the deep quench limit the small wavelength behaviour 
dominates and according to an earlier comment in section 
unci we expect to find /3 ^ TOq ^ for any dimension. For 
small values of m/mQ the asymptotic expressions of fd 
and gd allow us to calculate analytically the first order 
corrections of /? as a function of to/toq. In this way we 
find 



• d=l 



d = 2 



d = 3 



4 32 log 2to 
1 

mg TTTOq 



Too 



3 log 2 - 2 

log(TO/TOo) 



(34) 



(35) 



— (2tt/V3 - 7r(2 - n/V3)m/mo + ...) 
Too V / 

w Too 1(3.6276 - 0.584967TO/mo + ...) 



(36) 




the logarithmic corrections could render comparison with 
data difficult. In 3d the contribution of nonzero but small 
k modes causes a small shift of the numerical factor from 
4 to 27r/x/3« 3.6276. 



IV. MASSLESS Id THEORIES 



In section IIII Al we saw that for d = 1 and to = 
the propagator does not become stationary. However 
the situation is different when we consider a physical Id 
quantum system, for the following reason. A massless 
free theory is not physically meaningful. The infrared 
divergences impose the introduction of interaction coun- 
terterms of all orders in perturbation theory over the in- 
troduced coupling constant. The field renormalization 
finally results in the physical field defined as the exponen- 
tial of the original gaussian field (the vertex operator) . 
Thus in a physically meaningful Id system, interaction 
terms must always be present and the correlation func- 
tion is given by the expectation values of vertex operators 
^git}0(x)g-tt}0(x )^ Jqj. appropriate value of the constant 
q. This can be evaluated readily using the well-known 
property of gaussian integrals 



/^iq<f,{x) ^-iq<f,{x')\ 



(37) 

where C{x — x') ~ {(j){x)(j){x')) is the free propagator we 
have already found. From (|15p we obtain 



Uiq<t>{0,t) ^-iq(t>(r,t)\ 



-q^mot/i 



if r > 2i, 
if r < 2t. 



(38) 



m 



Thus the linearly increasing time dependence of C(r, t) 
leads to an exponentially decaying correlation function 
outside of the horizon and a static form inside the hori- 
zon. Therefore thermalization also occurs in Id systems. 
This has been shown to be the case for any massive to 
massless quench on a Id bosonic system, using the map- 
ping to the slab and the powerful methods of conformal 
field theory^i^. 



V. ANHARMONIC COUPLED OSCILLATORS 
(INTERACTING FIELD THEORY) IN 
SELF-CONSISTENT APPROXIMATION. 

Let us now consider a system of anharmonic coupled 
oscillators. The simplest form of a hamiltonian describ- 
ing such a system is 



FIG. 3: Effective temperature as a function of the final mass 
Pmo = Fd{m/mo) in units of the initial mass mo = 1. Inset: 
Asymptotic behaviour for small m. Notice the logarithmic 
corrections in 2d. 



In Id and 2d the k = momentum mode dominates 
so that the first order term is /3 = 4/too. In 2d however 



H 



(39) 



In the continuum limit this corresponds to the sim- 
plest form of an interacting quantum field theory, the 
(/)^ model. At t = we instantaneously change the mass 
from Too to to and at the same time the coupling con- 
stant from Ao to A. As before, we assume that initially 
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TABLE I: Asymptotic behaviour of the quench integral fd{s) 



d 


exact 




s ^ 1 


s — >■ CXJ 


1 

2 
3 


[2 log s + (Vl — s^/s) arccos s]/4 

[2(s - 1) - Vs^ - 1 arccos(l/s)]/4 

[(1 — s^)/2 — s^logs — s\/l — arccos s]/4 


(7r/2s + 21ogs)/4 

-(logs)/4 

(l-7rs)/8 


(.-1)76 
(s- 1)712 
(s- 1)712 


(logs)/4 

(l + 7r/4)s/2 

(log2-l/2)s74 



TABLE II: Asymptotic behaviour of the thermal integral gd{s) (|32|) 



d 


exact 


s ^ 


s — >■ oo 


1 




(7r/2s) + (logs)/2 




2 


-log(l-e-»)/s 


-(log s)/s 




3 




(7r76s7[l-(3./7r)] 


„-3 / 777 —3/2 

e n / As ' 



the system lies in the ground state of the hamiltonian 
before the quench. 

Such a model is non-integrable and can be solved only 
approximately. In this paper we will focus solely on the 
Hartree-Fock or self-consistent approximation. Roughly 
speaking in this approach we assume that the quar- 
tic interactions can be approximated by a 'mean field' 
quadratic term with a parameter that should be calcu- 
lated self-consistently. More specifically the 4>'^ inter- 
action term of the hamiltonian can be substituted as 
follows^! 

where we have taken into account that (0) = and the 
numerical factors are derived by Wick's theorem as the 
number of combinations of operator contractions. 

Such a substitution is justified in the large- A'^ limit of 
the linear cr-model, which is a variant of the 0^ model 
where the field has A^-components 

N 

^ = E E i'^' + + + iM^f^' (41) 

i—1 r 

In the limit — > 00 with AA^ kept fixed, the Hartree- 
Fock approximation becomes exact. 

Just by staring at (|40|) we notice that the second term 
corresponds to a mass term but with a 'mass' that has to 
be determined from the two point correlation function. 
The first term is just a number and does not affect the 
equations of motion but, as shown in appendix|Xl ensures 
the conservation of the total energy. Therefore we can 
define an effective mass mcff according to 

k 

and since the right hand side also depends on mgff, this is 
in fact a self-consistency equation for mes- Note that the 
effective mass should correspond to the pole of the cor- 
relation function on the imaginary axis in the complex 



fc-plane, which is what is physically measurable as the 
mass of the particles of the system. As a result of this ap- 
proximation, our initial non-integrable problem has been 
effectively reduced to an integrable and in fact free one, 
subject to the self-consistency equation. In our out-of- 
equilibrium case we should keep in mind that the effective 
mass will be time dependent. 

We will use two slightly different methods in applying 
this approach. The first one is a perturbativc method. 
After introducing the Schwinger-Keldysh method which 
is suitable for out of equilibrium problems, we soon real- 
ize that the usual perturbativc expansion does not con- 
verge and a resummation of Feynman diagrams using the 
Dyson equation is needed. This leads us to a simple 
ansatz for the asymptotic form of the two point correla- 
tion function at large times. The second method empha- 
sizes on the time evolution of the system and is based on 
a direct integration of the equations of motion in their, 
simplified by the self-consistent approximation, version. 
In order to solve these equations we employ an approx- 
imate analytical and an exact numerical method. The 
results of both calculations are in agreement with each 
other and additionally they verify our earlier ansatz. 

Before we start, it is worth to remind ourselves of the 
large- A^ results for the ground state of our system, since 
in order to proceed to the out-of-equilibrium problem we 
will need to know more about the initial properties of 
the system. This will also introduce us to a discussion of 
the renormalization procedure and its application to the 
present problem. 



A. Divergences and renormalization 

The initial two-point correlation function of our system 
in the large- A^ limit is simply that of a free system with 
a mass equal to its effective value 

"^e^ffo = -o + yE(<^') (43) 

k 
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The sum in the right hand side of (|43|) represents the 
fluctuations of the field. In the continuum limit this cor- 
responds to the integral 



(2n¥ 



1 



(44) 



Vffo 



which exhibits ultraviolet (UV) divergences in all dimen- 
sions. In Id and 2d these can be absorbed completely by 
a mass renormalization, while in 3d an additional cou- 
pling constant renormalization is required^. 

The mass renormalization amounts to allowing the 
bare mass mp to be divergent so as to compensate the 
divergent integral. The (finite) renormalized mass is de- 



fined by m-g/j. 



Stuq is 



SuIq where the mass countcrtcrm 



Orrin = — 



An f d'^k 



2^/WT^ 

The effective mass in terms of mo/j is then 



(45) 



'■offo 



d'^k 



Ao 

2 ./ (27:)^ 



2jk^ 



"cSO 



2v^ 



(46) 



which is finite in Id and 2d. 

In 3d there is still a logarithmic UV divergence in (|46|) 
which can be absorbed by a coupling constant renormal- 
ization. A suitable renormalization countcrterm can be 
determined by studying the 4-point correlation function 
and turns out to be of the form 



(5 An = 



1 



d^k 



(27r)3 8(P -f 777,2^)3/2 



(47) 



The resulting renormalized coupling constant Xqr satis- 
fies 



An 



A 



OR 



1 — Ann SXo 
and replacing in (|46p we obtain 



(48) 



Aoi? f d^k 



'offo — '"-OR 



2 J (2-)^ 



'offo 



'offo '"Ofl 



2v/FT^ 4(^+ 7772^)3/2^ 



(49) 



which is indeed finite. Note that in all dimensions the 
solution to the above equations is 



777offo = 777ni? 



(50) 



i.e. the renormalized mass is identical to the effective 
mass. In what follows we should keep in mind the well- 
known renormalization group result that in 3c? the critical 



point of this model corresponds to zero coupling constant 
i.e. in the continuum limit the macroscopic behaviour 
of the theory is effectively free. Therefore interactions 
are not physically meaningful in the continuum limit. In 
physical systems however, the existence of a finite lat- 
tice spacing that induces a natural UV cutoff renders all 
momentum integrals finite and there is not such a restric- 
tion. 

After the quench, the change in the mass and the cou- 
pling constant results in a change of the corresponding 
counterterms as well. This is required, otherwise new 
divergences in the equation for the effective mass are in- 
evitably born. On the other hand, this should not be 
regarded failure of renormalization theory as the 

latter does not have to apply to expectation values taken 
in states which are not obtained by rcnormaliscd oper- 
ators acting on the vacuum and our initial state is not 
such. In absence of a definite rule for the selection of the 
renormalization counterterms like the one for the ground 
state or thermal expectation values, several choices can 
be applied^^. In the present work we will be using the 
ground state counterterms of the theory after the quench. 

Since the field fiuctuations right after the quench are 
exactly the same as before it, the equation for the effec- 
tive mass right after the quench is 



A f d'^k 



1 



2 7 {2Tr)d2^k^+m, 
or introducing the mass renormalization 



OR 



(51) 



A f d'^k 



1 



1 



{2t:Y \2y^WT'. 



^OR 



2V^ 



(52) 



where 777^ is the renormalized mass after the quench and 



6m = — 



A f d^k 



2 J {2Tr)d 2 y/k^ + m 



(53) 



is the corresponding mass countcrterm. As before the last 
expression is convergent in Id and 2d, but not in 3d. If 
we use a coupling constant renormalization countcrterm 



SX 



1 



d^k 



(27r)3 8(/fc2 + to2j)3/2 



(54) 



we find 



2 /n+^ 2 .^R [ ( 1 

)="^-+Ty (2^ i^v^ 



'OH 



1 777,2^(0+) - 7772 

2^k^+ml 4(fc2 + 777^)3/2 



(55) 



which is only convergent in the trivial case 
™off(0^) — TT^oR where there is no jump in the ef- 
fective mass i.e. no quench at all. Recalling our previous 
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remark, we realise that this problem is due to the fact 
that the presence of interactions does not make sense 
in the continuum limit. In lattice systems however 
there is not such a problem and the practical meaning 
of the above is simply that TOeff(0~'') is very large. In 
the following we will therefore keep a large UV cutoff 
A in all expressions for the 3d case and investigate the 
dependence of our results on this. 

An interesting first observation is that as defined by 
([5^ the initial mass-square m^jj(0+) can be negative. 
Indeed for uiq < m, m^jj(0+) is always positive, but for 
mo > m the mass shift induced by the interactions is 
negative and if A is large enough then 771^^(0"'") < 0. In 
particular if m = the latter is always true. From a phys- 
ical point of view this negativity means that the quench 
can effectively drag the system into an unstable initial 
state like that of a double well (or generally 'mexican 
hat') potential. We will come back to this aspect of the 
problem later. 

The integration in ([5^ can be done analytically. Ex- 
pressing the integral in dimensionless form we have 

"^cff(0"^) + ^-^^rn^(r^hd{mR/mQ) (56) 

where is the total solid angle in d dimensions (fii = 
2,172 = Stt, = 47r) and the function hd{s) is defined 
as 

A— >oo 

hd{s) = I k'^-^dk ( , ^ , ^ ] (57) 



and can be easily shown to be 

r(iogs)/2 ifd = i, 

/id(s) = < (s - l)/2 ifd = 2, (58) 

[(s2 _ l)(logA)/4 ifd = 3. 

B. Perturbative approach 

In order to study the evolution of the effective mass for 
t > we have to calculate the two-point correlation func- 
tion and the usual way to do this is to use perturbation 
theory. The two-point correlation function is 

C(r,ti,t2) = («'o|0(O,ti)</)(r,i2)|*o) (59) 

where, as before, |4'o) is the ground state of the ini- 
tial hamiltonian. For simplicity let us first assume that 
Ao = 0, i.e. that there is no interaction before the quench 
so that |\E'o) is the ground state of a free hamiltonian. At 
this point we encounter an important difference with the 
usual QFT methods: at zero temperature the starting 
point of such a calculation is usually the following for- 
mula 

(0| r{0.(O,ti)0,('^,i2)exp (-1 /+^ dtH„t{t))}\0) 
(Q|r{exp(-z/+^dti7„,i(t))}|0) 

(60) 



But in our case this expression is inappropriate since it 
relies on the condition that |0) is the ground state of 
the free part of the hamiltonian and the interactions are 
swithed on and off adiabatically. In a quantum quench 
this is not valid because |5'o) is the ground state of a 
different hamiltonian and the changes are done instanta- 
neously. Thus we have to trace back to the origin of ([50)) 
which follows from the interaction picture formalism 

{'^o\T{(f>i{0, ti)4)i{r, 12) exp (— 

(61) 

where t is integrated over a contour JC that starts from 
some initial time i^, passes through ti and ^2 where the 
interaction picture field operators are placed, extends 
to some final time tf and then goes back to U so that 
times on the second half of the contour are considered 
to be later than those on the first half (Fig. |4]). This 



JC 

FIG. 4: The Schwinger-Keldysh contour for a quantum 
quench. 

is the well-known Schwinger-Keldysh method for non- 
equilibrium quantum systems^^— and is applicable to 
any choice of initial state. 

If l^'o) = |0) and the interaction is switched on and 
off adiabatically then we can extend ti —00 and 
tf — )• -|-oo. In this case, from the adiabatic theorem, 
the action of exp (+i dt Hint{t)) (i.e. the evolution 
operator along the second half of the contour) on |^'o) 
yields just a multiplicative constant and (j6ip reduces to 
(|60p. In the present problem we need to use the original 
expression (|6T|) instead. The initial time can be set to be 
ti = 0, when the interaction is switched on. However the 
same choice can be used even in the general case when 
Ao ^ 0, i.e. when the interaction is present before the 
quench, since as explained above, in our approximation 
the initial state is still that of a free theory but with the 
mass replaced by its effective value. 

It is worth to remark that an alternative way of deriv- 
ing the Keldysh contour is by using the slab construction 
mentioned earlier. In this approach one would have to in- 
tegrate in imaginary time from —L/2 to i.e. from 
one to the other boundary of the slab, then analytically 
continue the arguments of the operators from imaginary 
to real times as in r — >■ it and finally take the limit L — > 
thus recovering (|6ip . 

We can now expand (j61|) in powers of A. Accord- 
ing to the above, the zeroth order perturbative term 
(^o| 0i(O, ii)(?!)i(r, ^2) l^'o) is exactly the quench propaga- 
tor (jl3p with the masses mo and m replaced by their 
renormalized values. 

The first order correction C^^-'(ii,t2; k) corresponds to 
the single loop Feynman diagram (Fig. [5]). After applying 
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FIG. 5: First order Feynman diagram. 

Wick's theorem we find that C^^-'(ti,<2; k) reads 

^ J^dt'C{h,t';k)C{t2,t';k) J ^C{t\t';k') (62) 

The cxphcit form of the loop momentum integral 
/ d'^k C{t, t; k) is 



,d, /(t^ofc-'^fe)^ ml-rn^ 1 
d k I — cos2ajfei 



(63) 

but we also have to take into account the mass renormal- 
ization, which amounts to subtracting the UV divergent 
Feynman part and substituting the bare mass m by the 
renormalized mn 

d'^kC{t,t;k,m) ^ I d'^k (c{t,t]k,mR) ) 

J V '^^k,mRj 

(64) 

For brevity we redefine m, mg to be the renormalized 
masses mn, moR in all subsequent equations. Then 
can be written explicitly as 



d, I {^ok - Wfc)^ n^o - 



d'^k 



4a;^a;ofc 4iu^wofc 



cos 2uJkt ) (65) 



From the terms that remain in (j65p . the first one which 
is the time independent part is always convergent since 
it decays like k~^ for large k. On the other hand, the 
second term which is the time dependent part decays like 
k~^ cos 2uj}^t, which means that it converges in Id and 2d, 
while in id it is divergent only at t = 0. We also note 
that (p5)) does not suffer from infrared (IR) divergences 
in the massless case m = except in Id. 

Having analyzed the convergence of the loop integral, 
let us now calculate it. If we assume that m ^ 0, then 
the time-independent part has been calculated exactly 
already in section IIIIDI it is equal (up to a numerical 
factor involving the total solid angle in d dimensions) 
to mQ~^ fd{m/mo) where /(j(s) is given in Table HI On 
the other hand we recall that the time-dependent part 
has been shown to decrease with time. More specifically 
using the stationary phase method we find that for large 
times it decays like 



(to^ — ■m'^)m'^ ^ cos(2mi -|- (p) 



However for small times, this same time-dependent part 
can be important (or even divergent as we saw that hap- 
pens in Sd). 

Thus we are naturally led to the question whether it 
is safe or not to completely ignore the time-dependent 
part of the loop integral in calculating C^^^ for large 
times. If this is correct, then the effect of the loop dia- 
gram for large times is simply a mass shift equal to the 
time-independent part (recall that a mass renormaliza- 
tion counterterm induces a similar shift in the mass, but 
an infinite one). Higher orders in perturbation theory 
correspond to more loops and therefore one needs to em- 
ploy a resummation of all orders in order to compute the 
actual mass shift. Such a resummation can lead to a 
non-perturbative dependence of the mass shift and the 
correlation function on the coupling constant. Indeed if 
we calculate C^^^ from (|62]) assuming that the loop inte- 
gral is constant, then we find that the first order correc- 
tion increases linearly with time, i.e. it will eventually 
become larger than the zeroth order term and therefore 
the perturbative series does not converge. The required 
resummation can be done using the Dyson equation as 
described in the next section. 



1. Resummation using the Dyson equation 

As well-known the Dyson equation is an integral equa- 
tion satisfied by the two-point correlation function of an 
interacting theory that expresses the fact that the latter 
can be constructed from the propagator in a recursive 
fashion, using a number of 'skeleton' diagrams as building 
blocks. In our problem and in the mixed representation 
the Dyson equation can be written in the form 

G{ti,t2;k) - 

C{ti,t2]k)+ ( dt' j dt" C{ti,t';k)Y.{t' ,t"-k)C{t" ,t2]k) 

JK JK 

(67) 

where we denote the full correlation function by C and 
S is the self-energy insertion, i.e. a two-leg insertion also 
constructed recursively by the skeleton diagrams. 

In the large- TV limit the loop diagram is the only skele- 
ton diagram and therefore the self-energy is simply a loop 
of the full correlation function. The Dyson equation then 
takes the simplified form (Fig. [6]) 

C{ti,t2;k) = 

C{tiM;k)+ [ dt' C{ti,t';k)^{t')C{t',t2]k) (68) 

JK 



where 



(69) 



mo 



{mty/^ 



(66) 



taking into account the mass renormalization. Notice 
that comparing with (|42l) we realize that I](t) is nothing 
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but the shift in the mass-square 



'off 



(70) 



where C*(fc;mo,m) is the stationary part of the propa- 
gator. 

One can also check whether the remainder of the large 
time asymptotic form of the Dyson equation 




FIG. 6: Diagrammatic representation of the Dyson equation 
in the large- A*' limit. 

As we see the Dyson equation contains C explicitly in 
the right hand side but also implicitly in the definition 
of E. Thus it is difficult to solve in general. In many 
cases it is useful as a check of validity for an ansatz: we 
assume a particular form for C, substitute in the Dyson 
equation and check the consistency or determine any free 
parameters. This is how we are going to use it in our 
problem. 

Let us therefore construct an ansatz based on the hy- 
pothesis that the time-dependence of the loop be neg- 
ligible. Then the same can be assumed for the self- 
energy since this is nothing but a dressed loop, i.e. 
the sum of all 'cactus-diagrams'. This would mean 
that E(t) can be replaced by its large time station- 
ary value E* = limt_i.+oo E(t) or. according to ([70| . 
that the effective mass itself can be considered as time- 
independent and equal to its large time stationary value 
m* = limt_>.+co i^csit)- In other words we suppose that 
the effective mass simply jumps at the time of the quench 
from niQ to m* in which case the correlation function 
should simply be equal to the quench propagator for a 
quench from mo to m* . Note that our assumption is 
twofold: first we assume that rUeS tends to a station- 
ary value and second that this happens fast enough to 
approximate its evolution by a jump. 

According to the above, our ansatz is that the two 
point correlation function C(ti,t2]k) is approximately 
the same as the propagator itself but with m replaced 
by an asymptotic effective mass m* 

C{ti,t2;k;mo,m) ~ C{ti,t2;k]mo,m*) (71) 

We expect this relation to be asymptotically exact for 
large times, when any memory of the initial evolution of 
the effective mass will have been lost. 

Now that we have an ansatz for the correlation func- 
tion we can use the Dyson equation to check its validity 
and determine the value of the free parameter m* . By 
substituting (|7ip into ((68|) and (|69l) and replacing E(i) by 
E*, we find that the Dyson equation is satisfied exactly 
when m* satisfies the self-consistency equation 



E* = - 
2 



C*{k;mo,m*) - 

ZuJk 



dt' C{ti,t';k;mo,m){J:{t') - T.*)C{t' M]k]mo,ra*) 

Ik 

(73) 

with 

^(*')-^/(0(^^^'^'^^"^-™*)-i) (^4) 

tends to zero as supposed to. This is however a cum- 
bersome calculation and will not be presented. We will 
later show an alternative way to study the time evolution 
and verify our ansatz, but for the moment let us focus 
on the self-consistency equation ([72]) and investigate its 
solutions. 



2. Self- consistent calculation of the mass shift 
Written explicitly the self-consistency equation ((72|) is 



d'^k 



(75) 

where ujI = Vk'^ + m* ^ . 

Once again some comments about the id case are 
due as ([75| contains a logarithmically divergent integral. 
Therefore a UV cutoff A is assumed and the solutions 
771* will depend upon it. As can be verified however, the 
small A behaviour of m* is not affected by A. By the 
way the A-counterterm ([M)) would successfully remove 
the current divergence yielding the finite equation 



Afl 
2 



d^k 



(27r)3 
2ujk^l 



*2 2 
777 — 777 



(76) 



but according to the discussion in section IV Al this is 
supposed to be correct only for A — and therefore 
provides no more information than (|75p with a cutoff. 

Going back to the general case, if we make the mo- 
mentum integrals dimensionless then the self-consistency 
equation can be written as 



A nd 



2 {2ttY 



7770 



hd 



777 
771* 



(72) 



(77) 

where fd{s) and hd{s) are the previously defined func- 
tions dHl]) and (f57| . 

The above equations can be solved numerically or even 
analytically in several asymptotic limits like for A — > or 
771 — >■ 0. Fig. I7I8I and show plots of the solutions 77i* as 
a function of A for several values of tti in Id, 2d and Sd, 
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while Fig. [To] shows m* as a function of m for A ^ oo in 
Id and 2d. A first important remark is that for m ^ 
and small A the first order correction m* — m is linear in 
A, while for m = this is not true. Instead m* depends 
on A in a non-perturbative way in this case. The first 
order corrections in A for to = are summarized below: 



for d = 1 



m* = for all A 



(78) 



In fact it is more correct to talk about the limit 
TO — > 0, since to can never reach zero in Id. In this 
limit, TO* follows to to zero like 



TO 



TOo7r/2 



21og(TOo/TO) + 1 — IGtt^to^/Atoq 



for d = 2 



TO 



for d = 3 



TO 



Atoo 
"2^ 



log(TOo/A) 



47^^/2 



(79) 



(80) 



(81) 



independent of the cutoff. 



On the other hand for large A and to, to* increases like 
TO* ~ to^/2too in Id, to* ~ Am/Tr in 2d while in 3d the 
large A result is cutoff dependent. In addition, in 2d and 
for TO = and A — > oo we find to* — > 0.24954... toq. 

Of particular interest are the 2d results for to = 0. The 
fact that TO* 7^ means that from the critical evolution in 
the presence of interactions, there always emerges a finite 
effective mass which lets the system become stationary, 
in contrast to the free case. 



C. Time evolution 



We saw in section IV Al that the initial value of the ef- 
fective mass-square mlg{0~^) can be negative, while our 
ansatz suggests that its asymptotic final value is always 
positive. It is therefore worthwhile to investigate the time 
evolution of the effective mass in more detail. Although 
this can be done in the context of perturbation theory 
as in the previous section, an alternative and rather sim- 
pler way is by integrating the equations of motion for the 
field operator 0. Since the exact equations are nonlinear, 
even if we were able to solve them the solution would 
depend on the initial operators 0(0), 0(0) in a nonlinear 
way, thus preventing a direct application of the initial 
conditions ([7]) as done in section Hi Al Fortunately in the 
Hartree-Fock approximation this obstacle can be circum- 
vented since the 0^ interaction term of the hamiltonian is 
substituted by a quadratic 'mean field' term according to 
PO)) . As explained in section |Vl this substitution reduces 



the interacting into a free problem with a time-dependent 
effective mass given by 



A 



2uJk) 



(82) 



thus yielding a linear equation of motion. 

Even after this simplification however the problem is 
not trivial. In the following two sections we will first 
apply an approximate method that leads to an analyti- 
cal solution for small values of the coupling constant and 
later derive exact equations for the evolution of the cor- 
relation function which we will integrate numerically. 



1. Quasi- adiabatic self- consistent approximation 

A common approximation that could provide a com- 
pletely analytical treatment is the adiabatic approxima- 
tion which is based on the assumption that TOeff(i) varies 
slowly in comparison with the fast oscillations that char- 
acterize the solution^. This is not a reasonable assump- 
tion though, since it is the solution itself that determines 
the time dependence of TOeff(t). However as we show be- 
low, one can establish an alternative argument leading 
to the same approximate solution. The latter becomes 
equivalent to our earlier ansatz (|7ip for small A and pro- 
vides a first idea of the qualitative behaviour of the solu- 
tion. 

Since our problem is now free, it can once again be 
decomposed into a set of independent harmonic oscilla- 
tors. Of course the time-dependence of the frequency 
of each oscillator involves a summation over the whole 
set of them, but for the moment it suffices to consider 
a single quantum harmonic oscillator with an arbitrary 
time-dependent frequency uj{t). The hamiltonian is 



(83) 



The equation of motion for the field operator evolving 
under Lu{t) is 



4> + uj^{t)(f> = 



(84) 



If the frequency varies with time very slowly (adiabati- 
cally) then u/ui'^ <C 1 and as well-known the solution is 
given by 



(f>{t) =0(0) J 4^ cos ( f uj{t')dt' ) + 



;(i) 



^(0) 



^/uj{t)uj{Q) \Jo 



sin / Lu{t')dt' 



(85) 



A detailed derivation of the above equation in the quan- 
tum case can be found in appendix |Bl 

Although, as we said, the adiabaticity condition does 
not apply to our problem because uj{t) may exhibit os- 
cillations with the same frequency as the solution, the 
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FIG. 7; Solutions of the self-consistency equation (|75p in Id. The plots show the effective mass m* as a function of the coupling 
constant A for several values of m in units of mo = 1. Notice that as m — )■ the effective mass tends logarithmically to zero for 
all A. 
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FIG. 8: The same plot in 2d. Notice that, in contrast to the Id case, as m — >■ the effective mass tends to a non-zero value for 
aU A > 0. 
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FIG. 9: The same plot in 3d. The curves show a weak (but 
increasing for increasing A) dependence on the cutoff A. The 
dashed lines correspond to A = 10* while the full ones to 
A = 10^ 
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FIG. 10: Effective mass as a function of m for A — s> oo in Id 
(blue line) and 2d (red line) in units of mo = 1. The dashed 
straight lines are for reference. Notice that as m — >■ 0, m* — s> 
logarithmically in Id, while in 2d m* — 0.24954. 
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condition lo/lo^ ^ 1 is also valid when the amplitude 
of the frequency oscillations is sufficiently small in com- 
parison with the average value. This happens when the 
coupling constant A is sufficiently small so that from ([5^ 
Wcff(0 ~ ™- III this quasi-adiabatic approximation we 
can still use the last expression (|85|) as the solution to 
our problem. 

Having found the time evolution of <j) wc can use 
the initial conditions to derive the correlation function 
which is all we need in order to find TOoff(t). Re- 
call that from © we have (0(O)7r(O) + 7r(0)^(0)) and 
(02(0)) = l/2wo, (7r2(0)) = wo/2. By a direct calcula- 
tion 



HO) 



ALOQUj{t)Lj{Q) AujQLo{t)uj{0) 



2 jlu:k{t')dt' 



(86) 

Now going back to the interacting field theory model, 
we coclude that the equal time correlation function for 
each momentum mode {4>\(t)) is given by (|86p with ajfc(i) 
corresponding to the time-dependent effective mass ((82|) 



i.e. i^f.{t) = +mlff{t). Therefore the self-consistency 
equation for TOoff(i) is 



'-eff 



it) 



E 

k 
,2 



•"Ok 



4W0A;Wfe(t)Wfe(0) 



4wofcWfc (t)wfc(O) 



(87) 



This equation enables us to extract physical informa- 
tion about the evolution of the system through its only 
parameter mcs{t). A first observation is that mcff(i) de- 
pends on an average value over all previous times. The 
initial value of the effective mass m^jj(0+) seems to be 
crucial for the time evolution. If m^ff(0"*") > and A — )• 
then meff (i) exhibits weak oscillations and the adiabatic- 
ity condition is satisfied for all times. At large times 
the argument of the cos increases like 2ujkt where Uk is 
the time average of ojk(t). Therefore we can apply the 
stationary phase method to show that the oscillations 
decay in time and rriceit) indeed tends to a stationary 
value given by 



-oj, 



Ok 



\4wo/cW^Wfc(0) 



-L) 



(88) 



If however 771^5(0+) < 0, the small k modes exhibit, 
at least at short times, exponential instead of oscillat- 
ing evolution and the adiabaticity condition is no longer 
satisfied. The latter is also true in the marginal case 

™yo+) = o. 

Although is not the same as the corresponding 
equation of our ansatz ([75]). they are in perfect agree- 
ment for A — > where the quasi-adiabatic approximation 
is correct. In the next section we will see that it is pos- 
sible to construct a system of differential equations that 
describe the time evolution of meg (t) exactly, thus allow- 
ing us to investigate the large A regime. 



2. Exact time evolution equations and numerical solution 

Let us go back to the problem of a quantum harmonic 
oscillator with a time dependent frequency, described by 
the hamiltonian (|55| and the equations of motions ([M)) 
and start from scratch. Inspired by the adiabatic solution 
, we assume a solution of the form^^ 



0(0 



1 



cxp 



(89) 



where fl(t) is a suitable function that we wish to deter- 
mine. Substituting into ((84|) we find that ((89|) is the exact 
solution if il{t) satisfies the equation 




(90) 



By comparison with the constant frequency case we can 
find that the appropriate initial conditions for f2(i) are 



(91) 



Notice that if the derivatives of co are much smaller than 
w itself, we reproduce the quasi-adiabatic limit where 
r2(<) = uj{t) to first order. 

Taking into account the general initial conditions for 
0(0), 7r(0) we have 



'17(0) 



+7r(0) 



J^n{t')dA (92) 



from which, using once again the initial conditions ([7t. 
wc find that the equal time correlation function is 



1 



1 



2uj{0)luo 



+ 



2a;(0)wo 



2/>(i')dt' 



(93) 



In fact the only difference with ((86)) is that aj(i) has been 
substituted with r2(t). The overall result is that instead 
of (|84p one has to solve another differential equation (|90p . 
The advantage is that the former is an operator equation 
while the latter is an ordinary equation and it is easier to 
deal with real or complex-valued functions than opera- 
tors, especially since we will have to solve it numerically. 

In our interacting problem, the equal time correlation 
function for each momentum mode (0^(t)) will be given 
as before by (|93| where ftk{t) is also a function of k. 
Note that fifc(i) itself does not have to be of the form 
(k^ + M'^{t)y/-^ but for large k it is asymptotically equal 
to w/j(t), which ensures that nothing has changed as long 
as the convergence of the integral in ([5^ is concerned. 

The system of equations , ([50]) and (15^ completely 
determine the time evolution of the system. Although 
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very difficult to deal with analytically, it can be easily 
integrated numerically by discretizing the (k, t) space and 
iteratively applying the following loop: 

1. calculate flk{t) for each k from (I^Hl) . 

2. calculate {4)\{t)) for each k from (^5)) . 

3. calculate rn^f^it) from the self-consistency equation 

(1521), 

4. move one step forward in time t ^ t + dt. 

Fig. [TT] shows typical plots of the time evolution of 
the effective mass. For ni^ff{0'^) > we see that the 
latter exhibits decaying oscillations around an asymp- 
totic stationary value. We observe that this is the case 
not only for small values of A as we proved using the 
quasi-adiabatic approximation, but also for large ones. 
Moreover we find that even when TO^ff(0+) < in which 
case the quasi-adiabatic approximation fails, rri^git) in- 
creases quickly and soon becomes positive to follow an 
oscillating evolution similar to the previously described 
one. The reason is that the exponential growth of the 
momentum modes with k^ < —m^gij,) leads to a fast 
increase of TO^fj(t) that brings it to positive values, ceas- 
ing the exponential growth and leaving only oscillating 
modes^l"— . 

The asymptotic value m* as numerically estimated 
from the above method is systematically compared with 
that derived by our ansatz in the next section. It is 
remarkable that they arc in perfect agreement for all 
choices of values for the parameters we studied. 
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FIG. 11: Typical plots of the time evolution of the effective 
mass as obtained numerically both in Id. (a) The first plot 
corresponds to parameter values (mo,m. A) = (1,2,10) that 
yield a positive value for m^ff{0'^). The effective mass exhibits 
oscillations of decaying amplitude ~ t~^^^ about an asymp- 
totic value that is accurately predicted by our ansatz m* . (b) 
The second plot corresponds to (mo,m, A) = (1,0.5, 10) that 
yield a negative value for m^fj(O^). The initial exponential 
growth brings m^ff to positive values and as before nicff tends 
to the value m* found with our ansatz. The 2d and 3d cases 
are similar. 



3. Comparison of the quasi-adiabatic and numerical results 
with our ansatz 



Let us recall our earlier ansatz for the correlation func- 
tion C{k,t) stating that the latter is the same, at large 
times, as that for a free theory with m replaced by the 
final effective value m* = meff(i — >■ oo) which we find 
self-consistently, i.e. 



Cans {k, t) 



2ujt 



1 + 

^Ofe 



cos (2w)!t) 



(94) 



On the other hand the quasi-adiabatic approximation 
gives 



C'qa(A;, t) 



1 



2ujk{t) 



2wfc(0)cJofc 



cos {2jluk{t')dt' 



(95) 



while the exact evolution in the Hartree-Fock approxi- 
mation of the problem, presented in the last section, is 



Ccyi(k, t) 



1 



2Vlk{t) 



1 



(cjfc(O) - ugkY 
2ujk{0)ujok 



(96) 



The last two expressions differ only in that Vlk{t) is re- 
placed by u}k{t) in Cqii{k,t). An important difference 
between both last two expressions and Cans is that in the 
latter replaces uJk{0). Furthermore although the ar- 
gument of the cos in Cqa should tend to 2a;^t as in Cans 
for large <, this is not necessary for Cox- Thus Cans is not 
apparently consistent with cither Cox or Cqa, except for 
A — > where all of them are in agreement. 

Lacking an analytical argument to verify our ansatz, 
we rely on the numerical evaluation of the exact ex- 
pression and determination of the corresponding asymp- 
totic value of nies- Fig. [121 1131 show plots of the shift 
= ?n^ff(oo) — as a function of A for various choices 
of the parameter values and dimensionality, always in 
units Too = 1. The plots are based on the predictions 
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2 4 6 8 10 



FIG. 12: Comparison of numerical data (crosses) with our 
ansatz (lines) for m = 0. The plots are Eoo = m^ff(oo) — 
as a function of A in units mo = 1. The red line corresponds 
to 2d and the blue one to 3d with A = 100. 




FIG. 13: Comparison of numerical data (crosses) with our 
ansatz (solid lines) and quasi-adiabatic predictions (dashed 
lines) for several values of m (again in units mo = 1). The 
red lines correspond to 2d and m = 2, the green ones to 2d 
and m = 5 and the blue ones to Id and m = 2. It is clear 
that the numerics agree with our ansatz rather than the quasi- 
adiabatic approximation which is only good for small values 
of A. 

of our ansatz, of the quasi-adiabatic approximation and 
estimates drawn from numerical integration of the exact 
equations. By comparison we observe that the numerical 
data agree with our ansatz very well even for large values 
of A. In the contrary they do not agree with the quasi- 
adiabatic results, apart from first order in A. We conclude 
that, although our ansatz is not manifestly consistent in 
form with the exact solution, it however reproduces the 
exact results very successfully. 



VI. CONCLUSIONS 

We studied the problem of a quantum quench in which 
we simultaneously change the mass and the coupling con- 
stant of an interacting system. We restrict ourselves 
to the time dependent Hartree-Fock approximation and 



make the plausible hypothesis that for large times the 
two-point correlation function is the same as the propaga- 
tor but with a mass shift. We verify the self-consistency 
of our ansatz and derive the asymptotic effective mass as 
a function of m, mo and A which is shown to be correct 
by numerics. We point out that if mofr(i) approaches 
its final value moff(oo) sufficiently quickly then in the 
Hartree-Fock approximation the composite quench of the 
mass and the coupling constant is essentially nothing but 
a simple quench of the mass from mg directly to mcff(oo). 
In this case our ansatz would be justified and its generic 
success is probably an indication that such a fast 'relax- 
ation' process is indeed what happens. 

Our findings show that effective thermalization, one of 
the highlights of quantum quenches in free and Id confor- 
mal systems, is also possible in interacting systems such 
as the present model. Furthermore it is enhanced in some 
sense by the presence of interactions, since it occurs un- 
der more general conditions than in free systems (that 
is even in 2d massless systems). This is because of the 
shift in the effective mass of the system induced by the 
interactions. As this is their only effect in our approxima- 
tion, the effective temperature is still given by the same 
relation as in a free model but with m replaced by m*, 
thus depending on the coupling constant. In particular, 
the effective temperature is still momentum dependent 
as in the free case, but this should not be surprising: 
as explained in the introduction and the main text, in 
diagrammatic perturbation theory the Hartree-Fock ap- 
proximation amounts to keeping only 'cactus-diagrams', 
i.e. Feynman diagrams that can be constructed solely by 
loops, and ignores the effect of collisions between quasi- 
particles with different momenta that can induce a mix- 
ing of the different modes. The next order correction 
would be to take into account the 'sunset' diagram shown 
in Fig. [m 



k-q-q' 




FIG. 14: The 'sunset' diagram. 



We finally mention that, except for the stationary be- 
haviour, also the other qualitative features of the two- 
point correlation function that we observed in section 
nil Al for the mass quench in the free case, are general 
and present also in interacting models and for quenches 
of the interaction strength. This comment refers not only 
to the horizon effect for which it is obvious, but also to 
the characteristic oscillations^, either decaying or not, 
and is valid at least for integrable models (or even for 
sufficiently small deviations from integrability) since, as 
the present work suggests, the first effect to the quench 
is only a shift of the quasiparticle masses (equivalently of 
the poles of the scattering matrix) . 
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(|83|) . The latter can be diagonalized in terms of the 
instantaneous creation and annihilation operators a^{t) 
and a{t) defined by^ 



Appendix A: A conserved quantity 

In the Hartree-Fock approximation the effective fre- 
quency of each momentum mode is time dependent so 
that the time derivative of the corresponding single mode 
'hamiltonian' is not zero. Indeed 



dt 



hk{t) 



d_ 
'dt 



(Al) 



where in the last step we used the equations of motion 



0. 



However we can still construct a conserved quantity. 
From the self consistency equation ([5^ we sec that 



dt 



where 



Therefore 



dt 



hk{t) 



\C{t)4>l{t) 



(A2) 



(A3) 



(A4) 



and if wc take the expectation value on the initial states 
and sum over all momenta wc conclude that 



dh{t) - ^.Y^i^h.it)) = ]\C{t)C{t) = ^^(C2(t)) 



dt 



8 dt 



i.e. the following quantity 



h{t) ~ -XC'^{t) 



(A5) 
(A6) 



is conserved. As a demonstration of internal consistency, 
the last expression is precisely the Hartree-Fock form of 
the hamiltonian (p9|) according to the substitution (l40l) . 



Appendix B: The adiabatic approximation 

We consider the quantum harmonic oscillator with 
time dependent frequency, described by the hamiltonian 



and its hermitian conjugate. Notice that a{t) in the above 
relation depends on time only through Lu{t). The time 
evolution due to the dynamics of the problem is obtained 
from the Heisenberg equations of motion which in the 
case of operators that depend explicitly on time become 



da 



da 



-iLoa H a' 

2uj 



(B2) 



and its hermitian conjugate. The last equations form a 
system of linear differential equations that in matrix form 
looks like 



with solution 



Ait) 



2u 



(B3) 



t)-re.,{j:AitVt')(^^^^^ 



(B4) 



where T denotes time ordering. If the frequency varies 
only slowly (adiabatically) with time then uj/uj" <C 1 and 
A(t) can be approximated by 



Ait) 



-iu) 
+iuj 



(B5) 



which is diagonal, so that the solution to (jBSp is simply 



lit) = exp (-i /g w(t')dt')a(0) 



(B6) 



and its hermitian conjugate. Note that the first order 
correction due to the off-diagonal part of Ait) gives 

ait) ^e-'^o^'^^)dsaio)+ 
t 

/■^i'i:(!lg2»//^(.)d.^t(o) (B7) 

Zioit ) 



Keeping only the zcroth order term, we proceed to finding 
(/>(i) from 0(t) = (a(t) + a^t)) /^2a;(t) to obtau 
the main text. 
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